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ABSTRACT 

Using optical long baseline interferometry, we resolved for the first time the two wide 
components of HD 167971, a candidate hierarchical triple system known to efficiently 
accelerate particles. Our multi-epoch VLTI observations provide direct evidence for a 
gravitational link between the 08 supergiant and the close eclipsing O + O binary. 
The separation varies from 8 to 15 mas over the three- year baseline of our obser- 
vations, suggesting that the components evolve on a wide and very eccentric orbit 
(most probably e > 0.5). These results provide evidence that the wide orbit revealed 
by our study is not coplanar with the orbit of the inner eclipsing binary. From our 
measurements of the near-infrared luminosity ratio, we constrain the spectral classi- 
fication of the components in the close binary to be 06-07, and confirm that these 
stars are likely main-sequence objects. Our results are discussed in the context of the 
bright non-thermal radio emission already reported for this system, and we provide 
arguments in favour of a maximum radio emission coincident with periastron passage. 
HD 167971 turns out to be an efficient O-type particle accelerator that constitutes a 
valuable target for future high angular resolution radio imaging using VLBI facilities. 

Key words: Stars: early-type - Stars: binaries - Stars: individual: HD 167971 - 
Techniques: interferometry - Radiation mechanisms: non-thermal. 



1 INTRODUCTION 

Massive binaries made of O-type and/or Wolf-Rayet (WR) 
stars are known to be the scene of stellar wind collisions 
whose strength is intimately dependent on wind proper- 
tics. Beside t he hy drodynamic aspect of such a phenomenon 
(e.g. |Pittard]|2009l b colliding-winds offer the opportunity to 
study higher level p hysical processes such as copious ther- 
mal X-ray emission (IStevens et ai1ll992l ; iPittard k, Stevens] 
1 19971 : IPittard fc Parkin! I2010T)" o r even particle accelera- 
tion mechanisms (lEichler fc Usoylll993l; Bcnaglia fc Romerd 
120031 ; IPittard fc DoughertvlbOOfj ; iDe Beckerll2007r ). The lat- 
ter process considers a particle (electrons and hadrons) ac- 
celeration mechanism taking place in the presence of the 
strong shocks of the colliding-winds. Historically, the exis- 
tence of a population of relativistic electrons has been re- 
vealed by the detection of non-thermal (synchrotron) ra- 
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dio emission coming from several massive stars. More re- 
cently a significant non-thermal high energy counterpart 
has been confirmed in th e case of WR140 in hard X- 
ravs dSugawara et alj l201lT ) , and even in 7-rays for 77 Car 
|Farnier et alj|201l|) . 

So far, less than 40 massive stars revealed their ca- 
pability to accel e rate part i cles u p to relativistic energies 
(|De Beckerl 120071 : iBenaglial |2010T ). Even though it is now 
rather well-established that most of them are confirmed bi- 
naries, uncertainties remain concerning the multiplicity sta- 
tus of several targets. Among these, HD 167971 deserves a 
particular attention. This system is kn own to be one of the 
brightest synchrotro n radio emitters (jBieging et alj 1 19891 : 
iBlomme et al] | 2007fi. clearly establi shin g its particle accel- 
erator status. iLeitherer et alj (| 19871 ) and lDavidge fc Forbesl 
( 1988) showed that it consists of a close O + O eclipsing 
binary (P = 3.3213 d) with a third more luminous 081 star. 
Considering the strong opacity of stellar winds to radio pho- 
tons, the non-thermal radio emission is not expected to es- 
cape the wind-wind interactio n region within the close bi- 
nary (see lBlomme et al1 l2007. for a discussion), suggesting 
electrons should be accelerated in a putative interaction re- 
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Table 1. Observing log for HD 167971 and its calibrators. 



Epoch 


Date 


MJD 


Star 


Type 


# files 


Instrument 


Baseline 


1 


May 17, 2008 


54603.258 


HD 167971 


SCI 


11 


AMBER 


U1-U3-U4 


1 


May 17, 2008 


54603.275 


HD 179688 


CAL 


9 


AMBER 


U1-U3-U4 


1 


May 17, 2008 


54603.339 


HD 167971 


SCI 


9 


AMBER 


U1-U3-U4 


1 


May 17, 2008 


54603 356 


HD 179688 


CAL 


9 


AMBER 


TJ1-TJ3-TJ4 


2 


July 16, 2008 


54663.153 


HD 167971 


SCI 


8 


AMBER 


U1-U3-U4 


2 


July 16, 2008 


54663.170 


HD 179688 


CAL 


6 


A TV TTHITi 

AMBER 


TT1 TT(J TT 1 

U1-U3-U4 


2 


July 17, 2008 


54664.105 


HD 154088 


CAL 


5 


AMBER 


U1-U3-U4 


2 


July 17, 2008 


54664.120 


HD 167971 


SCI 


5 


AMBER 


U1-U3-U4 


2 


July 17, 2008 


54664.161 


HD 179688 


CAL 


8 


AMBER 


U1-U3-U4 


3 


March 29, 2011 


55549.389 


HD 167971 


SCI 


5 


AMBER 


A0-G1-I1 


3 


March 29, 2011 


55549.401 


HD 166521 


CAL 


5 


AMBER 


A0-G1-I1 


3 


March 29, 2011 


55549.417 


HD 167971 


SCI 


5 


AMBER 


A0-G1-I1 


4 


August 9, 2011 


55782.141 


HD 166521 


CAL 


5 


PIONIER 


A1-G1-I1-K0 


4 


August 9, 2011 


55782.151 


HD 167971 


SCI 


5 


PIONIER 


A1-G1-I1-K0 


4 


August 9, 2011 


55782.161 


HD 166191 


CAL 


5 


PIONIER 


A1-G1-I1-K0 



gion involving the third star. The study of Blomme et al. also 
revealed a radio light curve displaying long-term variations, 
compatible with a period of ~ 20 yr or more. HD 167971 ap- 
pears to be the brightest X-ray source in the NGC 6604 open 
cluster, with a rather hard thermal emission spectrum com- 
patible with pre-shock velocities to o large to be attribute d 
to the close binary colliding- winds (|De Becker et al]|2005h . 
With a period of about 3.3 d, the separation in the eclips- 
ing binary is too short to allow stellar winds to reach their 
terminal velocities before collision, leading therefore to low 
post-shock temperatures unable to produce a significantly 
hard X-ray spectrum. The latter point suggests again an in- 
teraction between the winds of the close binary with that 
of the distant 08 supergiant, although its gravitational link 
with the close binary has yet to be demonstrated. 

In this paper, we report on the results of long baseline 
interferometric observations of HD 167971 performed with 
the Very Large Telescope Interferometer (VLTI) at ESO 
(European Southern Observatory, Chile). The observations 
and data processing are briefly described in Sect .[2] The re- 
sults are discussed in detail in Sect .[3] and we conclude in 
Sect.H 

2 OBSERVATIONS AND DATA PROCESSING 
2.1 VLTI observations 

HD 167971 was a target of the VLTI on May 17, 2008 and 
on July 16, 2008, using 3 Uni t Telescopes (UT1- UT3-UT4) 
and the AMBER instrument |Petrov et al.l 120071 ). Our tar- 
get was also observed on March 29, 2011 with AMBER and 

3 Auxiliary Telescopes (ATs) in the A0-G1-K0 configura- 
tion, translating into baselines ranging between about 50 
and 120 m. The AMBER instrument was operated in the 
LR-HK mode, which provides a spectral resolution R ~ 35 
across the H and K bands. Calibrated visibility and closure 
phase values were computed using the latest public version 
of the AMBER da t a reduction pac kage (amd lib ve rsion 3.1, 
iTatulli et all 120071 ; IChelli. Utrera. fc Duvert]|2009h and the 



yorick interface provided by the Jean-Marie Mariotti Cen- 
tefl The data were processed fol lowing the same procedure 
as described by San a et al.l 1)20111 ). The accuracy of the abso- 
lute calibration was typically of 20% for the visibilities and 
five degrees for the closure phases. 

On August 9, 2011, HD 1 67971 was observed with the 
visitor-instrument PIONIER (|Le Bouquin et alj liolli ) fed 
by all four ATs of VLTI in the A1-G1-I1-K0 configura- 
tion. PIONIER was operated with a low spectral resolu- 
tion R « 15, that is, three spectral channels across the H 
band. Calibrated visibility and closure phase values were 
computed using the latest version of the PIONIER data re - 
duction package (pndrs version 2.3, Bouquin et al ]|201ll ). 
The accuracy was typically of 5% for the visibilities and one 
degree for the closure phases. Table [1] summarises all obser- 
vations obtained for the four observing runs, including the 
observations of calibrator stars used to estimate the instru- 
mental visibility and closure phase. 

Special care was taken to retrieve the absolute sign of 
the closure phases for both instruments, which were previ- 
ously calibrated with observations of known binaries. 

2.2 Binary modeling 

The angular separation between the two stars in the close 
eclipsing binary is expected to be of the order of 0.1 mas 
(assuming a circular orbit with a period of 3.3 d), and is 
therefore seen as an unresolved point source by our interfer- 
ometer. The 081 outer star is also seen as a point-like source, 
with an est imated angula r diam eter of 0.05 mas based on 
the work of iMartins et alj (|2005l ). We therefore a priori ex- 
pected to separate two point-like objects in our observations: 
the close binary and the 08 supergiant. For all four epochs, 
large closure phases and strong visibility variations with re- 
spect to wavelength, baselines and time are observed. This 
is illustrated in Fig.[T]for the closure phases at epoch 3. This 
suggests indeed a binary system with a flux ratio close to 

1 JMMC softwares: http://www.jmmc.fr 
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Figure 1. Illustration of the closure phases obtained during the 
epoch 3 observing run on AMBER (diamonds with error bars), 
and of the best-fit binary model (solid line). The x 2 value asso- 
ciated to this fit is about 0.8. 
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Figure 2. Reduced \ 2 map resulting from the fit of a binary 
model to the closure phases of the epoch 3 AMBER data set. 
The circle indicates the position of the sharply peaked minimum 
X 2 value. The figure appears in colour in the electronic version of 
the paper. 



unity. We therefore decided to model the system by two un- 
resolved stars with a constant flux ratio over the H and the 
K bands, a valid assumption for O+O binaries. 

We used the calibrated closure phases to derive the flux 
ratio and astrometry of the binary sys tem, following the 
procedure presented in lAbsil et al.l (l201ll ). In short, we con- 
structed binary models with various flux ratios and posi- 
tions within the interferometric field-of-view (~ 40 mas on 
the UTs, ~ 200 mas on the ATs at H band) , and computed 
the reduced x 2 goodness of fit for all models to populate a \ 2 
cube. The global minimum, its significance and its unique- 
ness were then determined by inspecting the cube. For all 
epochs, a unique binary solution was found. The best fit to 
the data for epoch 3 and the resulting \ 2 map are shown in 
Fig. [1] and Fig. [2] respectively. It is clear that the measure- 
ment of the binary separation is unambiguous: all locations 
are associated to \ 2 ^ 5 except the one identified with a 
circle in Fig.d which has X 2 ~ 0.8. 

The uncertainties on the best fit parameters were de- 
rived by inspecting the behavior of the reduced x 2 cube 
around the absolute minima (see lAbsil et al. I l201ll . for de- 
tails). The final error bars on the astrometry also take into 
account the estimated uncertainty on the spectral calibra- 
tion of the two instruments (i.e., central wavelength of the 
spectral channels). It was conservatively set to 3.2% in AM- 
BER (twic e the width of a spectr al channel) and to 2% in 
PIONIER (|Le Bouquin et al.l201ll ). Finally, we checked that 
a similar analysis, performed on the calibrated squared vis- 
ibilities instead of the closure phases, gives consistent re- 
su lts, although gene rally with larger error bars as discussed 
m lAbsil et al.l l|201ll ). The parameters of the best fit models 
for the four epochs are summarized in Table [2] 



3 RESULTS AND DISCUSSION 

3.1 First evidence for a wide orbit 

The main result of our investigation is a separation of 
HD 167971 into two components, respectively referred to as 
A and B. The brightest one, component A, should a pri- 
ori be associated to the 08 supergiant (see Sect . 13.21 for a 
detailed discussion). The best-fit relative positions of the 
two resolved components in HD 167971 are given in Ta- 
ble[21 an d plotted in Fig. [3] A second important result is 
the clear demonstration that the relative positions of com- 
ponents A and B changed significantly from one epoch to the 
other. In particular, it is noticeable that we detect a weak 
(but significant) change between observations 1 and 2, sep- 
arated by only about 2 months. These facts provide for the 
first time compelling evidence that the third star is gravita- 
tionally bound to the close binary, definitely confirming the 
triple status of HD 167971. Considering the lower limit on 
the orbital period su ggested by radio light curves (P ~ 20 yr, 
iBlomme et"ai]|2007l ), this strongly suggests a very large ec- 
centricity for the orbit with a periastron passage rather close 
to the epochs 1 and 2. 

The relation between periastron (d p ) and apastron (d a ) 
stellar separations, and eccentricity (e), is the following: 



On the basis of the values of the relative separation at epochs 
1 and 4, i.e., two extreme separations in our sampling, we 
derive an upper limit for the periastron over apastron sepa- 
ration ratio of about 0.6, translating into a lower limit on the 
eccentricity of about 0.25. It should be emphasized that this 
value is very conservative as the separation at epochs 3 and 4 
strongly underestimate the expected separation at apastron. 
The configurations at four epochs illustrated in Fig. [3] sug- 
gest an actual separation at apastron at least twice as large 
as at epoch 4, converting into an eccentricity of at least 0.5. 
Even though our analysis led to valuable constraints on the 
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Table 2. Relative position of component B with respect to component A, and near-infrared luminosity ratios for 
the three observations (average over H and K bands for AMBER, and over the H band for PIONIER). The angular 
separation (S) and the position angle (6>) are also given. 



Epoch 


Date 


<5RA 


S DEC 


S 


e 


Lb,hk/La,hk 






(mas) 


(mas) 


(mas) 


(dog) 




1 


May 17, 2008 


8.28 ±0.27 


-2.53 ±0.09 


8.66 ±0.28 


107.0 ± 1.1 


0.76 ±0.01 


2 


July 16-17, 2008 


8.78 ±0.28 


-1.03 ± 0.05 


8.84 ±0.28 


96.7 ±0.6 


0.81 ±0.04 


3 


March 29, 2011 


0.28 ±0.03 


15.03 ±0.48 


15.03 ± 0.48 


1.1 ±0.2 


0.83 ±0.01 


4 


August 9, 2011 


-1.53 ±0.03 


15.03 ±0.30 


15.11 ± 0.30 


354.2 ± 0.2 


0.75 ±0.01 



1 


1 1 1 1 1 1 1 1 
1 

I o component A i 


1 1 1 1 1 1 
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Figure 3. Relative positions of components A and B in 
HD 167971 at the four observation epochs. Component A, the 
brightest object, is arbitrarily located at coordinates (0,0). 



properties of the wide orbit in HD 167971, we note that the 
undetermined inclination, the uncomplete orbital coverage, 
and the lack of information from spectroscopic monitoring 
(especially close to periastron passage) inhibit our capabil- 
ity to derive the orbital parameters of the system. For this 
reason, we will refrain from speculating further on the or- 
bital parameters of the wide orbit and will come back to 
this issue a soon as the adequate observational information 
(interferometric and spectroscopic) will be made available. 

It is also interesting to note that our astrometric results 
for the wide orbit suggest a not so high inclination (defined 
as the angle between the line-of-sight and the normal to the 
orbital plane), although the eclipse measured for the close 
binary leads to an inclination angle for the short orbit of 
at least about 70° ([Davidee fc Forbes! Il98i ) . These results 
provide evidence for a significant inclination between the 
short and wide orbits in HD 167971, i.e., the orbits are not 
coplanar. This fact is important for the issue of star forma- 
tion srenaxi£s_dej2eiidJinTc^i initial cluster conditions (see 
e.g. ISterzik fc TokovinirJliooih . 



3.2 The components in HD 167971 

Our interferometric measurements alone could not lift the 
ambiguity in the association of components A and B to 
the close binary and to the distant supergiant. However, 
iLeitherer et alj (|l987l ) reported that the O supergiant is the 
most luminous and dominates the optical and UV spectra. 
Even though our measurements were performed in the near- 
infrared domain, the rather similar effective temperatures of 
the three stars in the system strongly suggest that the lumi- 
nosity ratio should go in the same direction in the infrared as 
in optical and UV domains. We will therefore consider that 
the supergiant is component A and that the close binary is 
component B. 

The luminosity ratio derived at every epoch can be dis- 
cussed in order to derive additional constraints on the spec- 
tral classification of the stars in HD 167971. Even though 
it is quite well-established that the third star should be 
an 08 supergiant, the spectral types of the components in 
the close binary are not fully determined. O n the basis of 
the fu ndamental parameter values given by iMartins et al.l 
(2005), we derived rough estimates of the expected lumi- 
nosity ratios assuming various spectral type combinations, 
and assuming the two components in the close binary to 
h ave the same spectral cl assification. According to the study 
of ILeitherer et a l. (1987), the members of the close binary 
should be main-sequence stars, with spectral types rang- 
ing between 05 and 08. The expected bolometric luminosi- 
ties and luminosity ratios, assuming an 081 third star (with 
Lbol ~4x 10 5 L©, IMartins et al.l [20051 . are quoted in Ta- 
ble|3] The values given in Table[2] for the measured lumi- 
nosity ratio seem to significantly change from one epoch to 
the other. It should be emphasized that the eclipsing nature 
of the close binary (component B) may be responsible for 
a significant change of its brightness as a function of time, 
mostly considering the depth and duration of the eclipses re- 
ported for instance by iDavidge fc Forbes! |l988). However, 
the propagation of the uncertainties on the ephemeris of 
the eclipsing binary over more than 2000 cycles (between 
the photometric and interferometric observations) prevent 
us from determining the orbital phase of the close binary 
at the epoch of our interferometric observations. In addi- 
tion, a slight difference may also appear between the three 
AMBER observations and the PIONIER measurement, as 
these instruments operate in different wavebands (HK and 
H, respectively). Assuming that luminosity ratios should be 
similar for both bolometric and K-band luminosities, our 
Lb,hk/La,hk ~ 0.8 measurements suggests spectral types 
for the stars in the close binary between 06 and 07 (typi- 
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Table 3. Potential stellar parameters for the stars in HD 167971, assuming main-sequence (upper part) and giant 
(lower part) luminosity classes for the components in the close binary. We assumed the close binary to be made of 
two identical stars. 



Sp. type 


05V 


05.5V 


06V 


06.5V 


07V 


07.5V 


08V 


Lbol (L©) 


3.2 x 10 5 


2.6 x 10 5 


2.0 x 10 5 


1.6 x 10 5 


1.3 x 10 5 


1.0 x 10 5 


0.8 x 10 5 


Lb /La 


1.60 


1.30 


1.00 


0.80 


0.65 


0.50 


0.40 


Sp. type 


05III 


05.5111 


06III 


06.5III 


07III 


07.5111 


08III 


Lbol (L©) 


5.0 x 10 5 


4.3 x 10 5 


3.6 x 10 5 


3.1 x 10 5 


2.7 x 10 5 


2.3 x 10 5 


2.0 x 10 5 


Lb /La 


2.5 


2.15 


1.80 


1.55 


1.35 


1.15 


1.0 



cally, 06.5 with an uncertainty of half a spectral type) for 
main-sequence stars, as illustrated in Fig. For a giant 
luminosity class, the close binary components should not 
be earlier than 08, assuming that the third star is at least 
as luminous as the close binary. On the basis of these es- 
timates we f avor the main-s e quenc e luminosity classes as 
suggested by iLeitherer et all (|l98?t ), and we consider that 
th e giant or even supergian t luminosity classes as claimed 
bv lDavidge fc Forbes! (|l98ct ) are unlikely. This result is also 
in agreement with the idea that the typical radius of evolved 
stars may be too large to fit within the small orbit of the 
eclipsing binary. C onsidering typ i cal ra dii and masses for 
06.5 stars given bv lMartins et al l (|2005h . and applying Ke- 
pler's law for a 3.3 d circular orbit, we estimate that super- 
giants would be too large and that about 80 % of the stellar 
separation would be occupied by the stars if they were gi- 
ants. If main-sequence stars are assumed, we estimate that 
only about 50 % of the stellar separation is filled by the stars. 

Assuming an (06.5V + 06.5V) + 081 spectral clas- 
sification for the stars in HD 167971, we estimate the to- 
tal visual absolute mag nitude My = 6.73 (f or typical visual 
magnitudes as given bv lMartins et al.ll2005h . For a V appar- 
ent m agnitude of about 7.35 as reported by ILeitherer et all 
(1987) between the eclipses, and a color excess E(B - V) of 
0.9, we derive a distance modulus of 11.29, translating into 
a distance of about 1.8 kpc. Considering the uncertainties, 
such a dist ance is in agreem ent with the distance of 1.7 kpc 
reported bv lReipurthl j2008h for the NGC 6604 open cluster, 
lending further support to the adopted spectral classifica- 
tion. 

We note that the additional information given here on 
the spectral classification of the components in HD 167971, 
with respect to previously published results, should be con- 
sidered with caution: stellar classification based on lumi- 
nosities are indeed generally less accurate than classifica- 
tion approaches based on spectroscopic criteria. This issue 
would therefore strongly benefit from a long-term spectro- 
scopic monitoring with application of state-of-the-art disen- 
tangling techniques to accurately classify individual compo- 
nents in the system, such as in the case of t he hierarchical 
triple system HD 150136 (|Mahv et al.l l2012h . provided the 
inclination of the wide orbit allows the measurement of sig- 

2 The uncertainty of 0.2 dex taken into account in this figure 
could seem overestimated considering the error bars given in Ta- 
ble[2]for the luminosity ratios, but we selected a comfortable con- 
fidence range in order to be sure to include the a priori unknown 
uncertainty on the theoretical luminosity ratio values, and there- 
fore get rid of any overintcrprctation of the results. 
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Figure 4. Expected luminosity ratios assuming various spectral 
types, along with main-sequence (V) and giant (III) luminosity 
classes for the close binary components. The long-dashed horizon- 
tal lines encompass the ratio values compatible with our interfer- 
ometric measurements, if one tolerates a 0.2 dex uncertainty. 

nificant radial velocity shifts of spectral lines of the three 
components. 

3.3 The radio emission in HD 167971 

The radio light curve presented by iBlomme et all (|2007h 
shows a bright non-thermal radio emission whose strong 
variability is worth discussing in relation with the long pe- 
riod eccentric orbit. The demonstration of a gravitational 
connection between the 08 supergiant and the eclipsing bi- 
nary provides a straightforwad interpretation scheme for the 
radio variability: the non-thermal radio photons should arise 
from an emission region coincident with the wind-wind inter- 
action region in the wide binary, where relativistic electrons 
(and particles in general ) are expected to be ac celerated. The 
light curve presented bv lBlomme et al.1 (2007) at 6 cm is the 
one that benefits from the best time coverage. Its maximum 
occurred close to 1988-1989, and its minimum between 1995 
and 1999. Assuming a period of about 20 years, the next 
maximum is expected to have occurred in 2008-2009, i.e., 
close to the time of our observations 1 and 2. In this context, 
it is tentative to attribute the maximum radio flux density 
to orbital phases close to periastron passage, i.e., when the 
separation is shorter. 

The non-thermal radio emission is indeed expected to be 
variable mainly for two reasons, respectively related to emis- 
sion and absorption. First, the intrinsic synchrotron mecha- 
nism responsible for the non-thermal radio emission is more 
efficient when the stellar separation is shorter (typically, 
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close to periastron passage) . Synchrotron emission is indeed 
dependent on the local magnetic field strength (B), where 
relativistic electrons are accelerated: L sy nch oc B 2 . Consid- 
ering the magnetic field strength depend ence on the dis- 
tance to the stars described for instance bv lUsov fc Melrose] 
(1992), B is expected to decrease significantly as the stel- 
lar separation increases, leading to a maximum synchrotron 
emission at periastron. In addition, one should keep in mind 
that the kinetic power injected in the wind-wind interaction 
region, and consequently the power available for particle ac- 
celeration, is expected to decrease as the separation increases 
(see discussion in Sect. 13. 4| ). Second, radio emission from the 
colliding-winds region is also expected to vary significantly 
because of the phase-dependent free-free absorption due to 
the stellar winds material along the line of sight. In very 
eccentric orbits, the absorption is expected to be stronger 
close to periastron, when the synchrotron emission region is 
more likely to be located within the r ~ 1 radiosphere. The 
combination of these two effect is especially striking in the 
case of the long period Wolf-Rayet + O-type binary WR 140 
l| White fc Becker! 1 19951 ). where the radio maximum occurs 
slightly before periastron, just before the strong increase of 
the free- free absorptiorjf]. In the case of HD 167971, how- 
ever, the impact of free-free absorption is a priori expected 
to be less spectacular: (i) the winds of the O-type stars in 
our target are not so efficient at absorbing radio photons as 
the WR wind in WR 140; (ii) HD 167971 is not so asymmet- 
ric as WR 140 and the stagnation point in the wide orbit is 
more likely to be located close to mid-separation. As a re- 
sult, we consider that the variability of the radio light curve 
in HD 167971 should be dominated by the variation in the 
synchrotron emission, with a maximum rather close to pe- 
riastron passage. 

Finally, considering the rather wide orbit of the sys- 
tem whose wind-wind interaction is responsible for the 
non-thermal radio emission, along with the large bright- 
ness of the synchrotron source in the system, we claim 
that HD 167971 would deserve high angular resolution radio 
imaging observations using Very Long Baseline Interferome- 
ter facilitites such as the VLB A, as already done in the cas e 
of a few systems such as WR 1 40 l|Doughertv et al" 2005h. 
WR1 46 (lO'Connor et al.1 l2005h. WR14 7 (|Williams et all 
Il997f). HD93129A dBenaglia et al.ll2010l) and CygOB2#5 
1 Ortiz-Leon et all 1201 ll ). This would allow to spatially dis- 
entangle the thermal and non-thermal contributions in 
HD 167971, and therefore more accurately quantify the syn- 
chrotron luminosity without contamination by the free-free 
(thermal) emission arising in individual stellar winds. 



3.4 Energy budget in HD 167971 

The results detailed in previous sections allow us to give a 
few words on energetic considerations related to the physics 
of massive binaries, and in particular to particle acceler- 
ation and non-thermal processes in general. The available 



energy for physical processes related to the wind-wind in- 
teraction comes from the kinetic power of the colliding 
stellar winds. For a stellar wind with mass loss rate M 
and terminal velocity 14o, the kinetic power is Pkin = 
0.5 MV^. For a long period orbit as the one revealed by 
our interferometric measurements, the stellar winds are ex- 
pected to have reached their terminal velocity before collid- 
ing. Considering the spectral types and luminosity classes 
discussed in Sect. 13.21 and assuming terminal velocities of 
30 00 km s~ , we used t he predicted mass loss rates given 
by iMuiires et al l d2012l) based on stellar parameters from 
iMartins et al. (2005) to compute kinetic powers of about 



1.9 x 10 36 and 0.7 x 10 36 ergs -1 respectively for components 
A and B. In order to evaluate the energy injected in the 
colliding- wind region, one has to consider the fractional solid 
angle of the wind-wind interaction region as seen from each 
star surface. Typically, one could consider that it amounts 
up to a few per cent. If a fraction of about 5 % is assumed, 
the total energy per unit time injected in the wind-wind in- 
teraction by both stars is about 10 35 ergs -1 . The efficiency 
of the conversion of collision mechanical energy to particle 
acceleration is expected to be of the order of one per cent, 
shared between electrons and hadrons (mostly protons), i.e., 
about 10 33 ergs -1 . The assumed value for the electron-to- 
hadron acceleration efficiency is still uncertain for colliding- 
wind binaries. Considering that this quantity is similar for 
environments accelerat ing cosmic-rays in general, it should 
be of the order of 1 % i|Longaidll992l ; IPittard fc Dougherty! 
2006). However, we should not reject the possibility that 
this efficiency ratio may be somewhat different, and possibly 
larger (1-10%). Consequently, at the level of approximation 
of this discussion, we will consider that 10 31 -10 32 ergs -1 
is potentially available for non-thermal emission processes 
involving relativistic electrons (i.e., so-called leptonic pro- 
cesses) . 

In massive star environments, the most efficient non- 
thermal leptonic process is expected to be inverse Compton 
(IC) scattering (i.e., energy transfer from relativistic elec- 
trons to UV and visible photons that are up-scattered in the 
high energy domain). Considering the ratio between mag- 
netic and radiation energy densities close to massive stars, 
IC scattering will be the dominant energy loss process for 
relati vistic electrons (jPittard fc Do ugherty 2006; [De Becker! 
12007? ) . Typically, one could expect that at most one per cent 
of the energy transferred to relativistic electrons is converted 
into synchrotron radiation in the radio domain. As a result, 

the synchrotron radio luminosity should be of the order of 
10 29_ 1Q 30 ergs -i 

This number is worth comparing to observation a l quan - 
tities as reported in the study by iBlomme et al. I (|2007h . 
To do so, we need radio flux density measurements at the 
same epoch in order to derive a spectral index. We selected 
in the data published by Blomme et al. measurements ob- 
tained in February 198qj, with measurements respectively of 
7.4 ± 0.4 mjy at 2 cm and 14.7 ± 0.3 mjy at 6 cm. From these 



3 We note that the situation may be somewhat more complex 
notably due to the Razin-Tsytovich effect , as illustrated in th e 
case of the massive binary CygOB2#8A llBlomme et al.ll2010l) , 
but the simple orbital phase dependence discussed here is correct 
to first order. 



4 That date was selected as it corresponds to the epoch closest 
to the maximum radio emission for which reliable flux density 
measurements exist at two wavelengths in order to determine a 
sepctral index. At the maximum (1988-1989), flux density mea- 
surements exist only at 6 cm. 
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quantities, we determine a spectral index a « —0.6, which 
is a comfortable non-thermal valu^fl We can therefore con- 
sider that the radio emission at the selected date is clearly 
dominated by synchrotron radiation, and that thermal con- 
tributions from the individual stellar winds contribute only 
a small fraction of the measured fl ux density at a given 
wavelength. Using equation (26) in |Pe Becken (|2007l ) and 
assuming a distance of 1. 7 kpc for the NG C6604 open clus- 
ter harboring HD 167971 (|Reipurth| [2008). we converted the 
flux density into an integrated radio luminosity of about 
10 30 ergs _1 . The latter value is close to the maximum ex- 
pected on the basis of the above discussion. This suggests 
that the electron-to-hadron particle acceleration efficiency 
is rather high in HD 167971, or that the energy injection 
in the particle acceleration mechanism is especially efficient. 
These considerations lend further support to the idea that 
HD 167971 is an especially valuable target for the investiga- 
tions of non-thermal processes in O-type multiple systems. 

These estimates suggest finally that the remaining en- 
ergy injected in the relativistic electron population (i.e., 
about 10 31 -10 32 ergs -1 ) should be available for IC scat- 
tering, likely responsible for a weak hard X-ray emission 
(and potentially soft 7-rays). Considering the sensitivity of 
current high energy facilities, the detection of such an IC 
scattering emission is unlikely, and next generation hard 
X-ray facilities su ch as ASTRO-H (J AX A mission to be 
launched in 2014, iMivata et al.l l201ll ) are required to de- 
tect and quantify this emission process, and therefore de- 
rive valuable constraints on the non-thermal and electron- 
to-hadron efficiency. 



4 CONCLUSIONS 

Our analysis provided for the first time evidence that 
HD 167971 is a gravitationally bound hierarchical triple sys- 
tem, with a very eccentric wide outer orbit. Using long 
baseline interferometry (with AMBER and P10N1ER in- 
struments at the VLT1), we resolved the system into two 
components: (i) the close eclipsing O-type binary and (ii) 
the additional O-type supergiant. The small fraction of the 
orbit covered by our observations suggests that the orbit 
between the two resolved components is wide, in agreement 
with the previously published lower limit on the orbital pe- 
riod of about 20 yr, on the basis of the radio light curve. Our 
results suggest an eccentricity of at least 0.25 (and most 
probably e > 0.5), even though additional interferometric 
observations are required to ascertain its value. 

Our measurements of the luminosity ratio allowed us to 
derive some constraints on the spectral classification of the 
members of the triple system. Considering that the brightest 
component is the 081 star identified by previous studies, our 
results favor an 06-07 main-sequence classification for the 
two stars in the close binary (assuming these two stars are 



5 A thermal radio emission component from a spherically sym- 
metric and homogeneous stellar wind is expected to be char- 
acterized by a spectral inde x close to 0.6, a s dem onstrated by 
IWright fc Barlow! lll975l) and lPanagia fc Fellil jl975h . A negative 
value for the spectral index is considered as a valuable tracer for 
non-thermal radio emission. 



very similar, in agreement with previously published photo- 
metric results). 

On the basis of these results, we can comfortably argue 
that the strong non-thermal radio emission associated to 
HD 167971 should come from the wind-wind interaction re- 
gion between the close binary and the additional supergiant. 
The rather wide orbit allows the bulk of the synchrotron 
radio emission to escape the stellar winds without undergo- 
ing strong free-free absorption. The large (though undeter- 
mined) eccentricity of the orbit provides a straightforward 
explanation for the strong variability of the non-thermal ra- 
dio emission in HD 167971. Assuming an orbital period of 
about 20 yr, as suggested by previously published radio light 
curves, our results provide hints for a maximum radio emis- 
sion occurring close to periastron passage. A more complete 
interferometric monitoring of HD 167971 is however needed 
to establish that coincidence with more accuracy. We em- 
phasize that our results lend additional support to the so 
called 'standard scenario' for non-thermal processes in mas- 
sive star environments, that considers efficient particle accel- 
eration in the colliding-winds region. Finally, we note that 
the wide orbit and the brightness of the synchrotron source 
in HD 167971 designate it as a potential good target for fu- 
ture high resolution radio imaging observations of colliding- 
wind massive binaries. 
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